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Abstract 

We evaluate the two-pion contribution to the nucleon electromagnetic form factors by use of dis- 
persion analysis and chiral perturbation theory. After subtraction of the rho-meson component, 
we calculate the distributions of charge and magnetization in coordinate space, which can be in- 
terpreted as the effects of the pion cloud. We find that the charge distribution of this pion cloud 
effect peaks at distances of about 0.3 fin. Furthermore, we calculate the contribution of the pion 
cloud to the isovector charges and radii of the nucleon. 



#1 email: hammer@itkp.uni-bonn.de 

#2 Address after Jan. 1 st , 2004: Institute for Nuclear Theory, University of Washington, Seattle, WA 98195-1550, USA. 
# 3 email : drechsel@kph.uni-mainz.de 
#4 email: meissner@itkp.uni-bonn.de 



1. That the pion plays an important role in the long-range structure of the nucleon is known since long. 
This can for example be deduced from the phenomenological analysis of the nucleon electromagnetic 
form factors (for an early attempt within meson theory see e.g. pQ). However, only with the advent of 
QCD and its spontaneously broken chiral symmetry, in which the pions emerge as pseudo-Goldstone 
bosons, this concept could be put on a firmer basis. Exploiting the chiral symmetry of QCD, the long- 
range low-momentum structure of the nucleon can be calculated within chiral perturbation theory 
(CHPT), which is the low-energy effective field theory of the Standard Model. For calculations of these 
form factors within CHPT, see [2]-|10|. Furthermore, since vector mesons play an important role in 
the electromagnetic structure of the nucleon, see e.g. ^^-^OJ , care must be taken in attributing a 
certain size or length scale to the various contributions (this is discussed in some detail in section 2 of 
|21j). A new twist to this picture was given in the recent paper of Ref. [22] > wnere an interpretation 
of the form factor data was given in terms of a phenomenological fit with an ansatz for the pion cloud 
based on the old idea that the proton can be thought of as virtual neutron-positively charged pion 
pair. A very long-range contribution to the charge distribution in the Breit-frame extending out to 
about 2 fm was found and attributed to the pion cloud. While naively the pion Compton wave length 
is of this size, these findings are indeed surprising if compared with the "pion cloud" contribution due 
to the two-pion contribution for the isovector spectral functions of the nucleon form factors, which can 
be obtained from unitarity or chiral perturbation theory. As it will turn out these latter contributions 
to the long-range part of the nucleon structure are much more confined in space and agree well with 
earlier (but less systematic) calculations based on chiral soliton models, see e.g. [2H]. Therefore it 
remains to be shown how to reconcile the findings of Ref. [22], based on a global fit to all nucleon form 
factors, with the results of dispersion analysis and chiral perturbation theory. 

2. First, we must collect some basic definitions. The nucleon electromagnetic form factors are defined 
by the nucleon matrix element of the quark electromagnetic current, 



{N(p')\qrQq\N(p)) = u{p') 



^ Fl (q 2 ) + ^-a^(p'-p) u F 2 (q 2 ) 
2m 



nip) , (1) 



with q 2 = (p 1 — p) 2 = t the invariant momentum transfer squared, Q the quark charge matrix, and m 
the nucleon mass. F\{q 2 ) and F2(q 2 ) are the Dirac and the Pauli form factors, respectively. Following 
the conventions of |17j . we decompose the form factors into isoscalar (S) and isovector (V") components, 



F l {q 2 )=Ft{q 2 )+r,F^{q 2 ) , i = l,2, (2) 

subject to the normalization 

if (0) = F?(0) = \ , F 2 S > V (0) = ^±^, (3) 

with K p (K n ) = 1.793 (—1.913) the anomalous magnetic moment of the proton (neutron). We will also 
use the Sachs form factors, 

G I E (q 2 )=F{(q 2 ) + ^Fl(q 2 ), G^q 2 ) = F{ (q 2 ) + F{ (q 2 ) , I = S, V . (4) 

These are commonly referred to as the electric and the magnetic nucleon form factors. The slope of 
the form factors at q 2 = can be expressed in terms of a nucleon radius 

,2,1 6 dFj(q 2 ) 

(r)i "if(0) dq 2 g2=0 ' (5) 
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and analogously for the Sachs form factors. The analysis of the nucleon electromagnetic form factors 
proceeds most directly through the spectral representation given 

Fi{q) = - / „ — o T~ ' l = 1 > 2 > I = S ,V, (6) 

in terms of the real spectral functions al{p 2 ) = ImF?(p 2 ), and the corresponding thresholds are 
given by /Iq = 3M T , /j,q = 2M n . Since the isovector spectral function is non-vanishing for smaller 
momentum transfer (starting at the two-pion cut) than the isoscalar one (starting at the three-pion 
cut), we will mostly consider the isovector spectral functions in what follows. We consider the nucleon 
form factors in the space-like region. In the Breit-frame (where no energy is transferred), any form 
factor F can be written as the Fourier-transform of a coordinate space density, 

F(q 2 )= / rf 3 re^p(r) , (7) 



with q the three-momentum transfer. In particular, comparison with Eq. © allows us to express the 
density p{r) in terms of the spectral function 

M = 4^ / . d » 2 — • (8) 



Mo 



Note that for the electric and the magnetic Sachs form factor, p(r) is nothing but the charge and 
the magnetization density, respectively. For the Dirac and Pauli form factors, Eq. Q should be 
considered as a formal definition. This equation expresses the density as a linear combination of 
Yukawa distributions, each of mass p. The lightest mass hadron is the pion, and from Eq. (JHJ) it 
is evident that pions are responsible for the long-range part of the electromagnetic structure of the 
nucleon. This contribution is commonly called the "pion cloud" of the nucleon and in fact this long- 
range low-g 2 contribution to the nucleon form factors can be directly calculated on the basis of chiral 
perturbation theory, as will be discussed later. 

3. Next, we wish to evaluate the two-pion contribution in a model-independent way and draw some 
conclusions on the spatial extent of the pion cloud from that. As pointed out long ago and further 
elaborated on |14j . unitarity allows us to determine the isovector spectral functions from threshold 
up to masses of about 1 GeV in terms of the pion charge form factor F n (t) and the P-wave ttttNN 
partial waves, see Fig. ^ We use here the form 

ImGUt) = -^7= |^(i)| 2 J+(t) , 
m\/t 

lmG v M (t) = -4= \F v (t)\ 2 J_(t) , (9) 



where q± = Wt/i — M 2 . The functions J±(t) are related to the t-channel P-wave 7rN partial waves 
f±(t) via f±(t) = F w (t) J±(t) in the conventional isospin decomposition, with the tabulated values 
of the Ji(t) from [21]. For the pion charge form factor F^ we use the standard Gounaris-Sakurai 
form ^5]. We stress that the representation of Eq. © gives the exact isovector spectral functions for 
AM 2 < t < 1QM 2 but in practice holds up to t ~ 50M 2 . It has two distinct features. First, as already 
pointed out in it contains the important contribution of the p meson (see Fig. ^) with its peak 

#5 We work here with unsubtracted dispersion relations. Since the normalizations of all the form factors are known, 
one could also work with once-subtracted dispersion relations. For the topic studied here, this is of no relevance. 
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Figure 1: Two-pion contribution to the isovector nucleon form factors. On the left side, the exact 
representation based on unitarity is shown, whereas the triangle diagram on the right side leads to the 
strong enhancement of the isovector spectral functions close to threshold. Also shown is the dominant 
p-meson contribution. 



at t ~ 30M 2 . Second, on the left shoulder of the p, the isovector spectral functions display a very 
pronounced enhancement close to the two-pion threshold. This is due to the logarithmic singularity 
on the second Riemann sheet located at t c = AM 2 — M^/m 2 = 3.98M 2 , very close to the threshold. 
This pole comes from the projection of the nucleon Born graphs, or in modern language, from the 
triangle diagram also depicted in Fig. ^ If one were to neglect this important unitarity correction, 
one would severely underestimate the nucleon isovector radii |15| . In fact, precisely the same effect is 
obtained at leading one-loop accuracy in chiral perturbation theory, as discussed first in (21 US] . This 
topic was further elaborated on in the framework of heavy baryon CHPT [3J |2j and in a covariant 
calculation based on infrared regularization [7j. Thus, the most important two-pion contribution to 
the nucleon form factors can be determined by using either unitarity or CHPT (in the latter case, of 
course, the p contribution is not included). 

We now want to separate the (uncorrelated) pion contribution from the p-contribution in the isovector 
spectral functions, that is we decompose the isovector spectral functions as 

Im Gj(t) =Im G^ 2n {t) + Im G V /{t) , I = E,M, (10) 

and analogously for Im F^fy)- Using Eq. (jHj), we can then calculate the pion cloud contribution to 
the charge and magnetization density in the Breit-frame. The p-contribution in Eq. (|10j) can be well 
represented by a Breit-Wigner form with a running width [Sj, 

Vn biMl^/tTJt) 

*° °r*m = (Mf^rm . '= E -**> (id 

with the mass M p = 769.3 MeV and the width T p (t) = g 2 (t - 4M^) 3 / 2 /(ASnvt), where the coupling 
g = 6.03 is determined from the empirical value T p (M 2 ) = 150.2 MeV, and the parameters bj can be 
adjusted to the height of the resonance peak. The corresponding expressions for the imaginary parts 
of the Dirac and Pauli form factors can be obtained from Eq. l|4"|). It is clear that the separation into 
the (uncorrelated) pion contribution and the p-contribution introduces some model-dependence. To 
get an idea about the theoretical error induced by this procedure, we perform the separation in three 
different ways: 
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Figure 2: The densities of charge and magnetization due to the pion cloud. Left panel: 47rr 2 / o(r) for 
the isovector Pauli (upper band) and Dirac (lower band) form factors. Right panel: A~Kr 2 p{r) for the 
isovector magnetic (upper band) and electric (lower band) Sachs form factors. 



(a) The two-pion contribution can be directly obtained from the two-loop chiral perturbation theory 
calculation of [Hj ■ Together with the ^-contribution of Eq. Qll[) , this calculation gives a very good 
description of the empirical spectral functions.* 6 We will use the analytical formulae given in [§] 
where the low-energy constant C4 was readjusted to avoid double counting of the /^contribution 
(see 123). 

(b) A lower bound on the two-pion contribution can be obtained by setting F n (t) = 1 in Eq. Q. 
This prescription does not only remove the p-pole but also some small uncorrelated two-pion 
contributions contained in the pion form factor. 

(c) To obtain the two-pion contribution, we can also subtract Eq. from the spectral function 
Eq. @ including the full pion form factor. # 7 The parameters bE = 1.512 and bu = 5.114 are 
determined such that the two-pion contribution at the p-resonance matches the two-loop chiral 
perturbation theory calculation of UJ. Variation of the bj around these values gives an additional 
error estimate. 

Using these three methods, we obtain a fairly good handle on the theoretical accuracy of the non- 
resonant two-pion contribution. 

4. We have now collected all pieces to work out the density distribution of the two-pion contribution 
to the nucleon electromagnetic form factors. Before showing the results, some remarks are in order. As 
stated above, the spectral functions are determined by unitarity (or chiral perturbation theory) only 
up to some maximum value of t, denoted t max in the following. Thus, we have simply set the spectral 
functions in the integral Eq. ((BJ) to zero for momentum transfers beyond the value tmax — 
In Fig. |2 we show the resulting densities for the isovector form factors weighted with Airr 2 . The 
contribution of the "pion cloud" to the total charge or magnetic moment is then simply obtained 
by integration over r. The bands reflect the theoretical uncertainty in the separation. For all form 

* 6 Note that on the right side of the p, the two-loop representation is slightly larger than the empirical one, so that we 
expect to obtain an upper bound by employing this procedure. 
* 7 A similar procedure was performed in [26) to extract scalar meson properties from the scalar pion form factor. 
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factors, the lower and upper bounds are given by methods (b) and (a), respectively. Method (c) 
generally yields a result between these bounds, except for the Dirac form factor where it gives the 
upper bound. The weighted densities for the isovector Dirac and Pauli form factors are shown in the 
left panel of Fig. [51 We see that these charge distributions show a pronounced peak around r ~ 0.3 fm, 
quite consistent with earlier determinations (see e.g. [231 128) ). and fall off smoothly with increasing 
distance. In the right panel of Fig. El we show the densities (again weighted with Anr 2 ) for the electric 
and magnetic Sachs form factors which come out very similar to the case of the Dirac and Pauli form 
factors. In comparison with Ref. [22], we generally obtain much smaller pion cloud effects at distances 
beyond 1 fm, e.g., by a factor 3 for p\{r) at r = 1.5 fm.# 8 

We have also studied the sensitivity of our results to the cut-off i max - While this may increase the 
value of the "pion cloud" contribution, it leaves the position of the maximum essentially unchanged. 
However, it is obvious from Eq. (jSJ that masses beyond 1 GeV and corresponding small-distance 
phenomena (r < 0.2 fm) should not be related to the pion cloud of the nucleon. Finally, we show the 
corresponding two-pion contribution to the charges and radii for the various nucleon form factors in 
Tabled The contribution of the pion cloud to the isovector electric (magnetic) charge is 20% (10%) 
in the model of Ref. j^2]. This is consistent with our range of values for the electric charge but a 
factor of 1.5 smaller than our lower bound for the magnetic one, see Tabled Furthermore, note that 
the pion cloud gives only a fraction of all form factors at zero momentum transfer. Normalized to the 
contribution of the pion cloud, the corresponding radii are of the order of 1 fm. In the model of |22j . 
these radii are considerably larger, of the order of 1.5 fm. Note that if one shifts all the strength of the 
corresponding spectral functions to threshold, one obtains an upper limit r max = y/3/2 M^ 1 ~ 1.7fm, 
assuming that the spectral functions are positive definite. 



*f(0) 


^(0) 


G V E (0) 




(r 2 )Y 


{r 2 )X 


(r 2 ) E 




0.07.. .0.08 


0.4... 1.0 


0.1. ..0.2 


0.4... 1.0 


0.1. ..0.2 


0.2... 0.3 


0.2... 0.3 


0.2...0.3 



Table 1: Two-pion contribution to charges and radii (in fm 2 ) for the various nucleon form factors. 
The radii are normalized to the physical charges and magnetic moments. 

5. In this letter, we have considered the two-pion contribution to the nucleon isovector form factors. 
The corresponding spectral functions can be obtained from unitarity or chiral perturbation theory 
from threshold to = 4M 2 - up to i max — 40M 2 in a model-independent way. Subtracting the contri- 
bution from the dominant p-meson pole, we have constructed the uncorrelated two-pion component 
(loosely spoken the dominant effects of the nucleons' pion cloud). To obtain an estimate about the 
theoretical error of such a procedure, we have performed this subtraction in three different ways. 
The corresponding charge densities obtained by Fourier transforming the spectral functions peak at 
distances of about 0.3 fm and show no structure at the larger distances. 
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* Note that our results are not in disagreement with the recent Jefferson Lab data on the ratio of the proton electro- 
magnetic form factors |29II8()| . The effect observed there occurs at momentum transfers beyond 1 GeV 2 and is thus not 
related to the pion cloud. 
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